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The effective Lagrangian of QED coupled to dyons is calculated. The re- 
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(N 
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f^ , suiting generalization of the Euler-Heisenberg Lagrangian contains non-linear 

^1* \ P and T noninvariant terms corresponding to the virtual pair creation of 

Qh! dyons. The corresponding P and T violating part of the matrix element for 

ix. 

Q , light-by-light scattering is considered. This effect induces an electric dipole 



moment for the electron, of order M~^, where M is the dyon mass. The 
current limit on the electric dipole moment of the electron yields the lower 
dyon mass bound M > 1 TeV. 
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I. INTRODUCTION 

Very precise measurements achieved during the last decade have opened up for a new 
approach in elementary particle physics. According to this, evidence of new particles can be 
extracted from indirect measurements of their virtual contribution to processes at energies 
which are too low for direct production. For example, the top quark mass as predicted from 
precision electroweak data |jl[| agrees to within 10% with direct experimental measurements 

i- 

This approach has recently been applied |^ for the estimation of possible virtual 

monopole contributions to observables at energies below the monopole mass. One-loop 

dyon- induced quantum corrections to the QED Lagrangian were discussed in [Q. Taking 

into account the violation of parity (and time-reversal symmetry) in a theory with monopoles 

0, the emergence of an electric dipole moment was first pointed out by Purcell and Ramsey 

0. More recently, the effect due to monopole loop contributions has been discussed [^|[. 

The calculation of quantum corrections due to the virtual pair creation of dyons is a 

very difficult problem because the standard diagram technique is not valid in this case. 

The difficulty is connected both to the large value of the magnetic charge of the dyon and 

the lack of a consistent local Lagrangian formulation of electrodynamics with two types of 

charge (see e.g. ^ and references therein). So, there is no possibility to use a perturbation 

expansion in a coupling constant. But one can apply the loop expansion which is just an 

expansion in powers of the Planck constant h. 

II. EFFECTIVE LAGRANGIAN 

It is known (see, e.g. [1^) that the one- loop quantum correction to the QED Lagrangian 
can be calculated without the use of perturbation methods. The correction is just the change 
in the vacuum energy in an external field. Let us review the simple case of weak constant 
parallel electric and magnetic fields E and H. We impose the conditions e|E|/m^ <^ 1 and 



e|H|/?7i^ <^ 1 such that the creation of particles is not possible. In this case the one-loop 
correction can be calculated by summing the one-particle modes — the solutions of the 



Dirac equation in the external electromagnetic field — over all quantum numbers |]IOl, [11 



For example, if there is just a magnetic field, H = (0, 0, H), the corresponding equation is 



[h'^id^ + ieAfj) - m]%l){x) 







(1) 



where the electromagnetic potential is A^ = (0, — ify, 0, 0). The solution to this equation 
gives the energy levels of an electron in a magnetic field fT^Jir 



w? + eH{2n - 1 + s) + fc^ 



(2) 



where n = 0, 1, 2 . . ., s = ±1, and k is the electron momentum along the field. In this case 



the correction to the Lagrangian is [10,12 



ALh = ^ I dk 
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[m 





~"ds 



2 + ^2)1/2 + 2 £(m2 + 2eHn + k^y/^ 

n=l 

[esH) coth (esH) - 1 - K^s^H^ 



(3) 



where the terms independent of the external field H are dropped and a standard renormal- 



ization of the electron charge has been made []T0[. 

If we consider simultaneously magnetic (H) and electric (E) homogeneous fields, then 
equation (0), as well as its classical analogue can be separated into two uncoupled equations. 



each in two variables [jT3|. Indeed, in this case we can take A^ = {Ez,—Hy,0,0) and the 
interactions of an electron with the fields E and H are determined independently. For such a 
configuration of electromagnetic fields the correction to the Lagrangian is (see [jlO|, P- 787) 

eH 



^L='£,j:Jdkelf\k). 

n=l n 



(4) 



Here elj^' is the correction to the energy of an electron in the combined external magnetic 
and electric fields, which is in the first order proportional to e^E'^. 

So, the total Lagrangian is L = Lq + AL, where Lq = (E — H )/2 is just the Lagrangian 
of the free electromagnetic field in the tree approximation, and can be written as 



oo \ _^2 tt2 

The logarithmic divergency can be removed by the standard renormahzation of the external 
fields and the electron charge: 

Ej-eg = -^3 E; Hrcg = Z^ H; e^eg = Z^, e, (6) 

where Z^ = 1 + (a/Svr) /(ds/s)e~™ * is the usual QED renormahzation factor. Thus the 



finite part of the correction to the Lagrangian AL' can be written in terms of physical 



quantities as (see |jTO[, P- 790) 

AL' = / ^e"""'' [{esE){esH) cot{esE)coth{esH) - 1] , (7) 



which in the limit E = reduces to the renormalized form of (|^). 

The series expansion of (|^) in terms of the parameters eE/rn^ <^ 1, eH/m? <^ 1 yields 



H' - E")' + 7(HE)' , (8) 



the well known Euler-Heisenberg correction [|T4 

3607r2m4 

where e^ = a. 

Let us consider how the situation changes if we consider the virtual pair creation of dyons 
in the external electromagnetic field. Using an analogy with the classical Lorentz force on a 
dyon of velocity v with electric {Q) and magnetic {g) charges ^ 

F = gE + ^H + vx (gH-(7E), (9) 

we shall assume that the wave equation for this particle in an external electromagnetic field 



can be expressed as [16, 



(z7^Z}^ - M)^(x) = 0, (10) 

where M is the dyon mass, and iD^ a generalized momentum operator, with D^ = d^ + 
iQA^ + igBf,. 



The potential A^ and its dual 5^ are defined by F^y = d^A^, — d^A^ = e^^pad'^B'' where 
Ffj_i, is the electromagnetic field strength tensoiQ and £0123 = 1- The potentials in the case 
of constant parallel electric and magnetic fields can be expressed as 

A>^ = {Ez, -Hy, 0, 0), B^ = {Hz, Ey, 0, 0). (11) 

It is easily seen that the solution to the equation of motion for a dyon in an external 
electromagnetic field can be obtained from the solution to the equation for an electron (|1]) 
by the following substitution 

eE^QE + gH; eH ^ QH - gE. (12) 

Using the same substitution as in Eqs. (§) and (|^, we obtain the following expression for 
the quantum correction to the Lagrangian, due to the vacuum polarization caused by dyons: 



12n^ J s 1271^ J s 





where a total derivative has been dropped. 

For the renormalization of this expression we can introduce the renormalization factors 

HI 



00 



'-1 _ 1 I V / "-^S ^-APs. ry~\ _ T y / "^^ ^-W^S 



12n J s ' ^ 1211 J s 





Z~' = l-^ -e-'''\ (14) 



which are generalizations of the definition Z^ of Eq. (|^). In this case the fields and charges 



are renormalized as [ 16,1 7 



-^reg — ^e ^g -^5 -^reg — ^e ^g ^ ] C^gg — Z^ZgC ] (yf^eg — ^e ^g 9 ■ (15) 

This relation (plSf) means that the vacuum of electrically charged particles shields the external 
electromagnetic field but the contribution from magnetically charged particles antishields 



it. This agrees with the results of [|T8| and [[19|] 



^ This definition is consistent only if OA^ = OB^ = 0, i.e., for constant electromagnetic fields or 
for free electromagnetic waves. 



Considering now the case of weak electromagnetic fields, the finite part of the Lagrangian 
AL', can, by analogy with (^, be written as 

AL' = ^^^AliQ' - gr + 7QV]i^' - ^r + [leg^/ + 7{Q' - /)^](he)^ 

+6QgiQ'-g')iEE)iE'-E')}. (16) 



The expressions (Jg) and ([Tq ) describe nonlinear corrections to the Maxwell equations 
which correspond to photon-photon interactions. The principal difference between the for- 
mula (^) and the standard Euler-Heisenberg effective Lagrangian consists in the appearance 
of P and T non-invariant terms proportional to (HE)(H — E ). It should however be noted 
that this term is invariant under charge conjugation C, since then both Q and g would change 
sign. 

If we consider separately the virtual creation of dyon pairs, then because of invariance 
of the model under a dual transformation (see, e.g. [§]) the physics is determined not by 



the values Q and g separately, but by the effective charge y/Q'^ + g'^. In the same way the 
operations of P and T inversions are modified. However we will consider simultaneously the 
contributions from vacuum polarization by electron-positron and dyon pairs. In this case it 
is not possible to reformulate the theory in terms of just one effective charge by means of 
a dual transformation. Moreover the Dirac charge quantization condition connects just the 
electric charge of the electron and the magnetic charge of a dyon: eg = n/2 whereas the 
electric charge Q is not quantized. 

It is widely believed, based both on experimental bounds and theoretical predictions 
pO| that the dyon mass would be large, M ^ m, where m is the electron mass. Thus, 
in the one-loop approximation the first non-linear correction to the QED Lagrangian from 
summing the contributions (|^) and (|16D can be written as 



where the P and T invariant terms corresponding to vacuum polarization by dyons have 
been dropped because they are suppressed by factors M~^. Thus, their contribution to 

6 



the effective Lagrangian will be of the same order as that of the ordinary QED multiloop 
amplitudes which we neglect. 

III. PHOTON-PHOTON SCATTERING 

Expression (0) yields the matrix element for low-energy photon-photon scattering. In 
order to determine it, we substitute into (|17D the expansion 

FA^) = (^ / d'? {l,Au - quA,) e''^\ (18) 

Corresponding to the second term of (|1^) we find 






ap 



(27r) 
where 



^ j d^qid'q2d%d''q,5{qi + ^2 + 53 + qi)A,{q{)A,{q2)A,{q:,)A,{q,) M^'^'^^ , (19) 



QgiQ'-g 



2^ 



^..p. = M>^^P^iq„ q„ q„ q,) = "^o'^^m^ ' ^^^^"^2 [qks " g'^'iqsq,)] • (20) 

Symmetrizing this pseudotensor one obtains the P and T violating part of the matrix element 
for light-by-light scattering. With all momenta flowing inwards, ki + k2 + k^ + k^ = 0, the 
matrix element takes the form 

M'liupa = l[M^,upaiki, h, h, h) + M^p^„{ki, k^^, /cs, k^j + M^^^p{ki, k^, /cg, ks) 
+M^pf,^{k2, ks, ki, k^) + My„^p{k2, k^, ki, k^) + Mp^^^^ks, k^, ki, ^2)] 

_ Q9[Q__2_9_J_( A«^LaL/3L(TLP I p MPLaLfTT/3,^ f^'' k'^ k-P k" k-f^ 

-e^r9'"'{hk,)ktk^2 - ej''g'"'{k2k,)k^,kl - e^/''g'"'{k2k,)k^k^, 
-ej''g^''{hk,)k^kl - ej"g^''{hk,)k'^kl - ej'' g^''{krk2)k'^kf' 



Since the interaction contains an e tensor, the coupling between two of the photons is 
different from that involving the other two, and the familiar pairwise equivalence of the six 



terms does not hold. The matrix element satisfies gauge invariance (with respect to any of 
the four photons), 

k'^M^.p^ih, k2, fcs, k^) = 0, etc. (22) 

We note that the above contribution to the matrix element is proportional to the fourth 
power of the inverse dyon mass, M' oc M~^. However, this result is only valid at low 
energies, where the photon momenta are small compared to M, being obtained from an 
effective, non-renormalizable theory. 

Thus, as a result of interference between two one-loop diagrams corresponding to loops 
with dyons and those with simply electrically charged particles there is an asymmetry be- 
tween the processes of photon splitting and photon coalescence [Q. The physical effect of 
this asymmetry will depend on the photon spectrum and the directions of the photon mo- 
menta with respect to the magnetic field. In particular, the asymmetry vanishes when these 
are perpendicular, i.e. for cos^ = 0. Furthermore, the asymmetry is linear in the product of 
the dyon charges, and proportional to the fourth power of the electron to dyon mass ratio. 

IV. ELECTRIC DIPOLE MOMENT 

The contribution of this matrix element (^Tj) breaks the P and T invariance of ordinary 
electrodynamics. Thus, among the sixth-order radiative corrections to the electron-photon 
vertex there are terms containing this photon-photon scattering subdiagram with a dyon 
loop contribution (see Fig. 1), that induce an electric dipole moment of the electron ||^^. 

Indeed, one can write the contribution of this diagram to the electron-photon vertex a^: 



This has been noted by LB. Khriplovich |21] — see also a recent paper by Flambaum and Murray 



^Of course, there are more diagrams. 



{p' + kiY — w? + ie {p — k'iY — fn'^ + ^e 
where Ma/3'yf^{ki, k2, k^, k) is the polarization pseudotensor representing the dyon box dia- 
gram contribution to the photon-photon scattering aniphtude, the low-energy limit of which 
is given by the pseudotensor M'^^^^ of Eq. (^) . 

In order to extract the electric dipole moment from the general expression (^) it is 
convenient, according to the approach by ||22| to exploit the identity 



Kf3^f.ih, k2, k3, k) = -r_M,^^,(A;i, k^, k^, k), (24) 

which can be obtained upon differentiating the gauge invariance condition of the polarization 
tensor [cf. Eq. (|2^)] with respect to k'^. 



Substituting (|2^ into (23) we can write the ee'f matrix element as 

M,,^{p\p,k) = e''{k)u{p')A^{p',p)u{p) 

= e^{k)k''uip')A,,{p',p)uip) (25) 

where e'^(fc) is the photon polarization vector and 

^.Ap^p) = -ji^sjd'k,d% ^ ^ ^ ^^M^,,Ah,h,k„k) 

^^ ^ / I / N2 ^~r~^ 7 — rv2 ^~r~'y ■ ^6) 

[p + /bi)^ — m^ + te [p — ks) — m'^ + le 



Since the matrix element (|25|) is already proportional to the external photon momentum /c, 
one can put fc = in A^,y after differentiation to obtain the static electric dipole moment. 



Then, following [^], we note that due to Lorentz covariance of A^j,, it can be written in 
the form 



KAp\ p) = [Ag^u + Ba^, + CP^-f, + DP,^^ + EP^P,) 75 + • • • (27) 

where we have omitted terms that do not violate parity, as well as those proportional to k^, 
and where a^u = {li^lu - ^vli^il'^, and Pf, = p^+ p'^. 

Substituting this expression into the matrix element Mee^yip'yP, k) of Eq. ( P5| ) one can see 
that there are two contributions to the P violating part, arising from the B and C terms. 

9 



In order to project out the dipole moment from (^), one has to compare Eq. p7| ) with the 
phenomenological expression for the electric dipole moment d^ [p^]: 

Meeyip^p^k) = e''{k)k''u{p')^^5cr^.u{p), (28) 

— * 

In the non-relativistic limit it corresponds to the interaction Hamiltonian —[de/2m)aE. 
Thus, multiplying ( p7|) by a^,y75 and taking the trace we have: 

rn 

4 = -^Tr[a^,75An- (29) 

In order to provide an estimate of the induced electric dipole moment we need to estimate 
A^'^. The first task is to evaluate the polarization pseudotensor Map-^^ corresponding to the 
virtual dyon one-loop subdiagram. If we were to substitute for M^p-y^ the low-energy form 
^'ap-fn of Eq. (|21|) into Eq. (]23|), we would obtain a quadratically divergent integral. 

On the other hand, straightforward application of the Feynman rules in QED with mag- 
netic charge (see, e.g., |T^) would give for the photon- by-photon scattering subdiagram in 
Fig. 10, 



M«/37m(^i, ^2, A^s, k:) = -— j d^q Tr (v. 



p- 



Here Fq represents the magnetic coupling of the photon to the dyon, which we take according 



to ref. IIT6| to be 



The vertex function depends on A;'', the photon momentum entering the vertex, and on n'^, 
a unit constant space-like vector corresponding to the Dirac singularity line. It was shown 



^It should be noted that the expression (|2l| ) contains contributions from such loop diagrams with 
all possible combinations of either three or one magnetic-coupling vertex Y p. 
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by Zwanziger p^ that although the matrix element depends on n, the cross section as well 



as other physical quantities are n independent. 

Calculations using this technique are very complicated and can only be done in a few 



simple situations |T6[, for example, in the case of the charge-monopole scattering problem 
P5| . We will here avoid this approach. 

While the integration over q in Eq. (|30|) is logarithmically divergent (the magnetic cou- 
plings in (|30|) are dimensionless) , after renormalization the sum of such contributions must 
in the low-energy limit reduce to the form given in Eq. (|2T|) . We also note that the sub- 
stitution of ( ^OD into Eq. (|26D yields a convergent integral. Thus, the following method for 
evaluating K^^ suggests itself. We divide the region of integration into two domains: (i) the 
momenta ki and k^ are small compared to M, and (ii) the momenta are of order M (or 
larger) . 

In the first region, the form (|21|) can be used, but since the integral is quadratically 
divergent, the integral will be proportional to M^. Together with the over-all factor M~^ 
this will give a contribution oc M~^. For large values of the photon momenta, the other 
form, Eq. (|3D|) can be used. This gives a convergent integral, and dimensional arguments 
determine the scale to be M~^. It means that 

The numerical coefficient has been estimated as l/(47r^)'^, one factor l/47r^ from each loop, 
and the 1/24 of Eq. ( ^9]) is assumed cancelled by a combinatorial factor from the number of 
diagrams involved. This is of course a very rough assessment. 

Now we can estimate the order of magnitude of the electron dipole moment generated 
by virtual dyons. It is obvious from Eqs. (p9|) and (0), that in order of magnitude one can 
write 

e'QgjQ' - g') m 
^^ {4^^ M^- (33^ 

This estimate can be used to obtain a new bound on the dyon mass. Indeed, recent 



experimental progress in the search for an electron electric dipole moment [^] gives a rather 
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strict upper limit: de < 9 ■ 10~^^e cm. If we suppose that Q r^ e, then from ( p3D one can 
obtain M > 2 ■ 10^ m ^ 10^ GeV. This estimate shows that the dyon mass belongs at least 
to the electroweak scale. 

The above estimate coincides with the bound obtained by De Riijula for monopoles, 
from an analysis of LEP data, but it is weaker than the result given in [^ , where the limit 
M > 10^ GeV was obtained. The authors of ref. used the hypothesis that a radial 
magnetic field could be induced due to virtual dyon pairs. In order to estimate the effect, 
they used the well-known formula for the Ueling correction to the electrostatic potential, 
simply replacing the electron charge and mass with those of the monopole. But the Ueling 
term is just a correction to the scalar Coulomb potential due to vacuum polarization and 
cannot itself be considered as a source of a radial magnetic field. Indeed, there is only one 
second order term in the effective Lagrangian that can violate the P and T invariance of 
the theory, namely AL' oc EH. But in the framework of QED there is no reason to consider 
such a correction because it is just a total derivative. The reference to the ^-term, used in 
to estimate the electric charge of the dyon, is only relevant in the context of a non-trivial 
topology (e.g., in the 't Hooft-Polyakov monopole model) where their limit applies. In this 
case there are arguments in favour of stronger limits on the monopole (dyon) mass (see, e.g.. 



One should note that the dyon loop diagram considered above can also contribute to the 
neutron electric dipole moment. The experimental value dn < 1-1 ■ 10~^^e cm p8[ will in 
the naive quark model with m ^ 10 MeV allow us to obtain an estimate of the dyon lower 
mass bound which is similar to the one obtained for the electron. 

Acknowledgements 

We acknowledge numerous useful conversations with A. Gazizov, V. Kiselev and G. 
Calucci. One of us (Ya. S.) is very indebted, for fruitful discussions, to Prof. I. B. Khriplovich, 
to whom belongs the idea of the above described mechanism of an electric dipole moment 
generated in QED with a magnetic charge. 

12 



Ya. S. also acknowledges support by the Alexander von Humboldt Foundation, and, in 
the first stage of this work, by the Fundamental Research Foundation of Belarus, grant No 
F-094. This research (P. O.) has been supported by the Research Council of Norway. 



13 



[1] 

[2] 

[3] 
[4] 

[5 

[6 

[7 
[8 
[9 



[10 



[11 



REFERENCES 

J. Ellis, G. L. Fogli and E. Lisi, Pliys. Lett. B333 (1994) 118. 

CDF Collaboration, F. Abe et al., Pliys. Rev. Lett. 74 (1995) 2676; 
DO Collaboration, S. Abaclii et al., Pliys. Rev. Lett. 74 (1995) 2632. 

A. De Riijula, Nucl. Phys. B435 (1995) 257. 

S. G. Kovalevich, P. Osland, Ya. M. Slinir and E. A. Tolkachev, Preprint ICTP 
IC/188/95, Trieste, 1995; ( |hep-th/960lT33| ). 



M. Sachs, Ann. of Phys. 6 (1959) 244; 

L. M. Tomil'chik, Sov. J. of Exp. and Theor. Phys. 44 (1963) 160. 

E. M. Purcell and N. F. Ramsey, Phys. Rev. 78 (1950) 807; 
N. F. Ramsey, Phys. Rev. 109 (1958) 225. 



V. V. Flambaum and D. W. Murray, |atom-ph/9604005 



P. Osland and Ya. M. Shnir, Bergen Univ. preprint, 1996-04, ( hep-ph/9606298 ) . 



V.I. Strazhev and L.M. Tomil'chik, Electrodynamics with a Magnetic Charge (Nauka i 

Tekhnika (in Russian), Minsk, 1975); 

V.I. Strazhev and L.M. Tomil'chik, Sov. J. Part. Nucl. 4 (1973) 78. 

A.I. Akhiezer and V.B. Berestetskii, Quantum Electrodynamics (Interscience Publish- 
ers, New York, London, Sydney, 1965). 

G. Wentzel, Quantum Theory of Fields (Interscience, New York, 1949). 



[12] V.B. Berestetskii, E.M. Lifshitz and L.P. Pitaevskii, Quantum Electrodynamics (Perg- 
amon Press, Oxford, 1982). 

[13] V.G. Bagrov and D.M. Gintman, Exact Solutions of Relativistic Wave Equations 
(Kluwer Academic Publishers, Dordrecht-Boston-London, 1990) pp. 103. 



14 



[14] W. Heisenberg and H. Euler, Z. Phys. 98 (1936) 714; 

V. Weisskopf, Mat. Fys. Medd. Dan. Vid. Selsk. XIV (1936) 6. 

[15] D. Zwanziger, Phys. Rev. 137 (1965) B647; 
S. Weinberg, Phys. Rev. 138 (1965) B988; 
A.S. Goldhaber, Phys. Rev. 140 (1965) B1407. 

[16] M. Blagojevic and P. Senjanovic, Phys. Rep. 157 (1988) 233. 

[17] G. Calucci and R. Jengo, Nucl. Phys. 223 (1983) 501. 

[18] C. Goebel and M. Thomaz, Phys. Rev. D30 (1984) 823. 

[19] E.A. Tolkachev and Ya.M. Shnir, Sov. J. Nucl. Phys. 55 (1992) 1596. 

[20] P. Goddard and D.Ohve, Rep. Progr. Phys. 41 (1978) 1361. 

[21] I. B. Khriplovich, private communication. 

[22] J. Aldins, S. J. Brodsky, A. J. Dufner and T. Kinoshita, Phys. Rev. Dl (1970) 2378. 

[23] I. B. Khriplovich, Parity Nonconservation in Atomic Phenomena^ (Gordon and Breach, 
1991). 

[24] D. Zwanziger, Phys. Rev. D6 (1972) 458. 

[25] A. Rabl, Phys. Rev. 179 (1969) 1363. 

[26] J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics, (McGraw-Hill, New 
York, 1964), p.l70. 

[27] J. P. Jacobs et al., Phys. Rev. A52 (1995) 3521; 
E. D. Commins et al., Phys. Rev. A50 (1994) 2960; 
W. Bernreuther and M. Suzuki, Rev. Mod. Phys. 63 (1991) 313. 

[28] I. S. Aharev et al., Phys. Lett. B276 (1992) 242. 



15 



FIGURES 




FIG. 1. Typical three-loop vertex diagram. The closed line represents a dyon loop. 
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